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A HOST-PARASITE PROBLEM 



INTRODUCTION 

Various biological processes in plants and animals occur as periodic func- 
tions of time, in which the rhythm of these occurrences is governed by a 
biological ' 'clock' 9 (1). Galston and Dalberg (2), for example, recognized 
that an inducible enzyme system could show rhythmic properties. In the 
simplest case, an enzyme is induced in the presence of a substrate which 
itself is destroyed by the induced enzyme. The substrate is produced at a 
constant rate, but the increasing amount of enzyme reduces the total sub- 
strate concentration. 

This substrate reduction, in turn, decreases the rate of induced enzyme 
formation, thus allowing the concentration of substrate to increase. Hence, 
a periodic behavior results between the amount of enzyme induced and the 
concentration of the substrate. 



An almost analogous condition exists in the case of the classical 
parasite* ' problem presented in this study. 



'hOSt- 



Printed in U.S.A. 464 



1 



© Electronic Associates, Inc. 1964 
All Rights Reserved 
Bulletin No. ALAC 64066 



PROBLEM DESCRIPTION 

A typical ' 'host-parasite* 9 situation can be des- 
cribed as follows: in a parasitic relationship in 
which the birth of a parasite (i.e., an increase in 
parasite population) causes the demise of its host 
(i.e., a decrease in host population), the rate of 
host population growth or decay depends on the 
number of living hosts and the number of encount- 
ers between host and parasite. 

The following assumptions were made: 

1. In the absence of parasites, the net rate of 
increase in host population is exponential, 
the food supply being unlimited. 

2. In the absence of hosts, the rate of parasite 
decrease is also exponential. 

3. The rate of encounter between hosts and 
parasites is constant. 

SYSTEM EQUATIONS 

The differential equations describing this host- 
parasite relationship are 



— = K.H - K 4 HP 



IT = " K 2 P + K 3 HP 



(1) 
(2) 



where H = host population as a function of time 
P = parasite population as a function of time 
Ki = overall growth rate of hosts 
K2 = overall decay rate of parasites 
Kg = K4 = rate of host parasite encounters 
t = time 

PROBLEM ANALYSIS 



Suppose, also, that the hosts increase at a rate of 
5% per hour (in the absence of parasites), and that 
the decay rate for parasites is 10% per hour (in 
the absence of hosts). Further, assume that there 
is one encounter between host and parasite per 
5000 hours. 



Therefore, 



K ± = 0.05, /hour 
K 2 = 0.10, /hour 



K 3 = 2 x 10-4, /host-hour 
K4 = 2 x 10" 4 , /parasite- hour 

COMPUTER PROGRAM 

The voltage- scaled computer equations are 
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and 
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Time Scaling: The potentiometer settings 

(Kx), (K 2 ), (1000K 4 ), and (2000K 3 ) in Equa- 
tions (3) and (4) are reasonably low. This sug- 
gests that the problem should be speeded up by a 
factor of, say, 10 to 1. Accordingly, a time scale 
factor of /3 = 0.1 is introduced to the denominator 
of all potentiometers going into integrators. 
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(6) 



The specific problem is to analyze the rhythmic 
behavior resulting between host and parasite in 
the situation described, through solution of the 
system differential equations. 



Equations (5) and (6) are the scaled computer 
equations from which the system model is con- 
structed. 



Solution of these equations on an analog computer 
requires numerical values. Suppose that, initially, 
there are 100 hosts and 200 parasites, and that the 
pattern of the relationship between them will be 
established by the time the host population reaches 
2000 and/or the parasite population reaches 1000. 



Please note that in addition to the 10:1 time 
scaling employed in this problem, there is also 
an implied time scaling which relates computer 
time (in seconds) to problem time (in hours). 
Thus, 1 second of computer time is equivalent 
to 10 hours of real time. 
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The computer diagram for this problem is shown 
in Figure 1. Tables showing the potentiometer and 
amplifier assignments are given in Figures 2 and 3. 
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Figure 3. Amplifier Assignment Sheet 



RESULTS AND CONCLUSIONS 

Figure 4 shows a plot of typical results obtained 
for a problem of this type. The rhythm of in- 
crease and decrease in both host and parasite 
populations, once established, continues as long 
as the same constraints are imposed. As shown, 
when the number of parasites is near zero (during 
the first 50 hours, for example), the number of 
hosts shows a rapid increase. This increased 
availability of hosts, however, provides an in- 
creased availability of parasite breeding ground 
and, during the second 50 hours, the parasite 
population makes an even more rapid increase to 
a peak occurring when the host population itself 
is near zero. 



Figore 1. Computer Diagram - - Host-Parasite Problem 
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Figure 2. Potentiometer Assignment Sheet 



Figure 4 
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In addition to this typical solution of the stated 
problem, a number of other interesting experi- 
ments can be performed with such an analog 
computer model. For example, if in a particular 
culture medium, it is known that ambient tempera- 
ture variations have an effect on either growth or 
decay rates, these coefficients can be varied 
simply through the manipulation of hand- set poten- 
tiometers. 

Another possibility for experimentation, at the 
expense only of a slight increase in the complexity 
of the equations, would be the introduction of a 
second type of parasite with different characteris- 
tics. 

It can be seen that use of an analog computer 
model in the solution of this problem gives the 
advantages of rapid solution, immediate availability 
of graphic results, and quick and easy parameter 
variations. 



, EQUIPMENT COMPLEMENT 

The major pieces of equipment required for this 
simulation include: 

• 2 integrators 

• 1 high gain amplifier 

• 2 inverters 

• 1 summing amplifier 

• 6 potentiometers 

• 1 quarter- square multiplier 

• 1 X- Y plotter. 
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